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Abstract 

A method of calculations of meson mass spectra in the framework of the model 
of quasi-independent quarks is discussed. Meson mass spectra evaluated with the 
. help of the Dirac equation with the quasi-Coulombic and confinement potentials, as 

I well as with the help of phenomenological mass formulae, are presented. Parameters 

' of the quasi-independent quarks model, which are calculated on the basis of mass 

values of pseudoscalar and vector mesons, are given. Possible relations between 
^ , parameters of the quasi-independent quarks model and constants of the Standard 

' model of strong and electroweak interactions are considered. 
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0> . 1. Introduction 

<^ , 

! It is well known that the calculation of hadron mass spectra on the level of experimental 

^ I data precision pQ still remains among the unsolved QCD problems due to some technical 

O I and conceptual difficulties, which are related mainly with the nonperturbative effects, 

^ • such as confinement and spontaneous chiral symmetry breaking. For instance, although 

recently a considerable progress in lattice QCD was achieved, at present the calculations 
Oh. of hadron mass values on the basis of the first principles of QCD demand too much 

r-| ! computation time especially for the light hadrons. So data interpretations and calculations 

of hadron characteristics are frequently carried out with the help of phenomenological 
models. Moreover the identification of hadrons bearing constituents, distinguishable from 
standard ones, now is the topical problem for strong interaction physics. For this purpose 
a detailed description of characteristics of standard hadrons is needed. Especially it 
concerns the evaluation of masses for ground and exited hadron states with light quarks 
and/or antiquarks. 

One of the most interesting and simple hadron models is the relativistic model for 
quasi-independent quarks. The main statement of this model is that hadron's properties 
can be described considering the hadron as a system of independent constituents (or 
quasi- independent ones with weak residual interactions), which move in some mean self- 
consistent field [2]. For the description of constituent's motion the one-particle equation 
(Schrodinger, Dirac or Klein-Gordon-Fock equation) can be used. The confinement of 
color particle is available by using of a linear rising potential, for which the hypothesis 
of its fiavour independence has been proposed [3l HI |5] . The relativistic model for quasi- 
independent quarks has been applied with the particular potential entered in the Dirac 
equation for the description of meson properties in Refs. P [7|. The potential used is 
a generalization of Cornell potential and consist of a relativistic vector quasi-Coulomb 
potential and a relativistic scalar linear rising potential. In the framework of this model 
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the hypothesis of universahty of confinement potential offered in Refs. [3], HI |5j has been 
confirmed for heavy and hght quarks. The coefficient of hnear rising scalar potential 
(string tension) has been found to be a = 0.20 ± O.OlGeV^^. Below we present the results 
of calculations of mass spectra of mesons and model parameters, which may be related 
directly to constants of Standard Model. In order to obtain the consistent treatment of 
1 and ^ mesons it is nessasary to take into account spin-spin interaction between 
quark and antiquark as well as the values of the meson mean field energy [7]. Note that 
the account of these terms are most important in the case of the pseudoscalar mesons. 
We made a link of quasi-independent quark model with the constituent quark model 
and propose that the values of the meson mean field energy belong equidistant energy 
levels. We present also mass spectra of excited meson states with 1 = 1 consisting of 
u—, (i— quarks and antiquarks with help of analytical mass formulae, which have the 
structure peculiar for mass formulae of the independent quark model. These formulae 
may be considered as the generalisations of the Chew-Frautschi and Nambu-Veneziano 
relations for the Regge trajectories of superlight mesons (mesons which consist from u—, 
d— quarks and antiquarks only). We compare results obtained with existing data and 
make predictions for mass values of unobserved meson states. Besides that, we consider 
possible structures of some mesons which are under discussion at present, such as scalar 
mesons and vector excited mesons. 

The report is organized as follows. We briefly review the basic equations and the main 
statements of the model in Section 2. Here the Dirac equation in the external potential 
with the Lorentz scalar and vector parts is transformed into a form suitable for numerical 
calculations. The values of model parameters and meson masses are presented in Section 
3. In Section 4 the phenomenological mass formulae for isovector radially and orbitally 
exited gg'— mesons are presented and the the results of calculations for mass spectra of 
exited mesons are displayed. In Section 5 the model of quasi- independent quarks is related 
to the constituent quarks model. We discuss the results obtained and the problems of 
identification of some meson states in vector and scalar channels in the last section. 

2. Basic statements of the relativistic independent quark model 



According to the main statement of the independent-quark model a hadron is consid- 
ered as a system composed of a few non-interacting with each other directly, valency con- 
stituents (quarks, diquarks and constituent gluons) having the coordinates r^, z = 1, A^, 
and moving in some mean field. We suppose that this field is a colour singlet confining 
field which is produced by the constituents and nonperturbative QCD vacuum and takes 
into account the effects of creation and annihilation of a sea of qq pairs as well. Further- 
more, to simplify calculations it is assumed that this mean field is spherically symmetric 
and is a quasi-classical object possessing some energy Eq inside a hadron. Each of 
constituents interacting with the spherical mean field gets the state with a definite value 
of its energy Ei, so that the hadron mass can be evaluated as 

Mh = Eo + E, + ... + En. (1) 

The wave function tpi for any constituent is a solution of a single-particle equation with 
the mean field static potential f/(rj). Angular dependence of the single-particle wave 
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functions in a stationary state for the spherically symmetric potential may be separated 
in a well-known manner, and in order to evaluate the constituent energies Ei it is necessary 
to solve the radial equations with the model potentials for each constituent. Note that it 
is impossible to evaluate Eq without additional assumptions, and this quantity take as a 
phenomenological parameter in the framework of the model. 

According to these assumptions the mass Mm of the qq - meson in the main approxi- 
mation can be evaluated as 



Mm = Eom + Ei{n[,ji) + ^sK, J2). (2) 

where Ei{n'[,ji) , i = 1,2, are the energy spectral functions or the mass terms for 
the i—th quark (antiquark) and represents the relativistic effective energy of the i—th 
quark/ ant iquark moving in the mean field inside the meson. Here n[ and ji are the radial 
quantum number and the quantum number of the angular moment correspondingly for 
the z— th constituent. 

The term Eom contains the contribution from the mean field energy and the possible 
nonpotential corrections, which cannot be evaluated in the frame of mean field approxi- 
mation. This term has a nonzero value only for some meson ground states. The terms 
Eiinl^ji), i = 1,2, which represent the energies of the constituents in the mean field, 
should be evaluated for quarks from the solution of the Dirac equation (for diquarks and 
constituent gluons from the solutions of Klein- Gordon- Fock equation): 

^/A~T^V^i(ri) = [{cxipi) + Piim, + V^) + Vi] V^,(ri), (3) 

with Ei{nl,ji) = ^/Xi + mf, i = 1,2, Vo{r) = ar/2 and Vi(r) = —2as/3r, where the 
model parameters a and as have meanings of the string tension and the strong coupling 
constant at small distances, correspondingly. As it is seen from Eq.Q, the addition of 
some constant to the scalar linear confining potential Vo(r) is equivalent to the addition of 
the same constant to the quark mass and vice versa, while the addition of some constant 
to the vector quasi- Coulombic potential Vi(r) is equivalent to the energy shift of opposite 
sign. 

It is well-known that the solutions of Eq. ([3]) with the total angular momentum j and 
its projection m can be represented as 



where n = r/r, the subscript i here and below is omitted. If A; = —uj{j + 1/2), where 00 
is an eigenvalue of the space-parity operator, the system of the radial Dirac equations for 
the fermion in the mean field with the definite energy sign and spin projection reads 

{V^T^ -Vo-Vi~ m)f = -^-^—^g - g', 

r 

{VX + m^ + Vo-Vi + m)g = ii±^/ + (4) 
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Using Eqs.(jll) one can derive the second order equation for the "large" component 
/(r), and then making a substitution 

r I 1 "^/^ 

(y9(r) = r/(r) Vb(r) — '\/i(r) + m + v A + 

one comes on to the model radial equation for y9(r) in the following form: 

k{k + 1) 3(\/o' - V[f 



+ 



4(VA 
(K," - V^i") 



+ m2 - \/i + \/o + m) 2( VaT^ - 1^1 + K) + 



(5) 



The ground states of standard gg'— mesons consist of the S-wave 1 and 0^^ mesons 
and the most extensive set of precise data [Ij for heavy mesons is related to the radially 
exited S-wave 1 mesons. For this type of mesons k = \ and the parameter A entering 
Eq.(l5]) can be calculated with the help of the radial equation: 
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or 



The right-hand side of Eq.([6]) has a singularity at the origin, and when r — > it 
behaves as 3/4r2 — 4a;^/9r2. Therefore one should keep < 3/2 in order to prevent a fall 
at the origin. At the infinity r — cxd the effective potential behaves as the oscillatory one 
and tends to a'^r'^ /A. The form (jH]) of the effective potential partly clarify some varieties 
of the confining potential used in potential models. We see, that for light quarks the 
confining potential is mainly oscillatory, while for heavy quarks it is linearly rising in the 
main approximation. 



3. Evaluations of model parameters and meson masses 



The model equation ([6]) can be solved only by numerical methods. For calculating its 
eigenvalues the computer code algorithm, which was based on the Numerov three point 
recurrent relation |8j|, has been used jointly with the regularization procedure for the 
singular potential at r = 0, which is analogous to the procedure worked out, for instance, 
in Ref . [9] . The accuracy criterion for fitted parameters is the maximum value of the errors 
for evaluated hadron masses, which should not exceed typical experimental errors 30 -r- 40 
MeV. Thus, the parameter errors written below are our estimation of systematic errors 
of the model, which have been found in this manner using mesons mass spectra data [1] . 

The brief description of parameter fitting is as follows. When calculating the values of 
the model parameters we do not suppose a priori the validity of the hypothesis of fiavor 
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independence for the confining potential. First of all the model parameters m, as and a 
for the bb and cc radial excitation of 1 states were found. The fits were carried out for 
each meson family independently. The existing experimental data allowed to find these 
parameters as well as to prove the validity of the model. In this manner it was found 
that the value of a is the same for the bb and cc states within the systematic errors of the 
model and equal to a = (0.20 ± Om)GeV^. 

Taking into account that there are no well established data for the higher radial exci- 
tations of 1 light mesons, the a value obtained was used for calculating m and as for 
light quarks. Note that when doing the calculations we neglect the isotopic mass splitting 
of the u- and d-quarks because it is beyond of the model accuracy. It should be noted that 
the terms Ei{nl,ji), i = 1,2, and i?om in Eq-(l2]) represent the main contributions to the 
energy of the system due to the interaction of the constituents with the mean field and 
the mean field energy, respectively. However, it is obvious, that in the spin independent 
potential the constituent energies Ei{nl, 1/2) are the same for the 1 and 0~^ mesons 
with the identical radial quantum numbers and the mass differences can be the result of 
the spin-spin interaction and the different values of the i?om terms. 

In order to estimate the spin-spin interaction between quark and antiquark we take 
into account the following expression for the interaction of this type which has been used 
in the relativized quark model [TO] : 

where Si and S2 are the spin operators and Ei and E2 are the energies of quark and 
antiquark, respectively. Evaluation of the expectation value of this operator in the first 
order of perturbation theory for the wave functions of quark and antiquark shows that the 
spin-spin interaction depends mainly on the Ei and E2, the total spin of quark- antiquark 
system and the strong coupling constant as- So one can use the modified mass formula 
for S-wave 1 and 0^+ mesons, which takes into account spin-spin interaction between 
quark and aniquark, in the following form [7]: 

= Eom + Ei{n[, 1/2) + 1/2) + 4 < S1S2 >q,^, Vss (8) 

Although in some cases for vector heavy mesons the contributions of spin-spin inter- 
action are less than the model precision, nevertheless, for pseudoscalar mesons for all 
bound states of quark and antiquark, the spin-spin contributions should be taken into 
account. The magnitude of the spin-spin interaction constant Vss is of the order of 100 
MeV for light quarks and of the order of 10 MeV for heavy quarks, while the values of the 
Eom parameter are: Eo{K) = -llSMeV, Eo{7]c) = -230MeV, Eo{J/<i/) = -125Me1/, 
Eo{7r) = -236Mev, Eo{rib) = -QOOMeV, Eo{T) = -AbOMeV, Eo{Ds) = -lUMeV, 
Eo{Bc) = —270MeV . The values of i?om for mesons, which were not pointed out above, 
within the model accuracy may be set to zero. 

We obtain the following values for the coupling as and the quark masses: 

m = (0.007 ± 0.005) GeV, = 0.65 ± 0.15, 
= (0.14 ±0.03) GeV, = 0.47 ± 0.10, 

rric = (1.28 ± 0.05) GeV, = 0.33 ± 0.03, 

b 



mb = (4.60 ±0.10) GeV, a° = 0.27 ± 0.02. 
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Within the model accuracy all results for the well known experimental data on the 
1 and 0"+ mesons, which are composed of the quark and antiquark with u-, d-, s-, c- 
or 6-flavours, are consistent with the presented above model parameters. The evaluated 
mass spectra of the pseudoscalar and vector mesons are displayed in the Table 1. Thus in 
the framework of the model the flavor independence of the confinement potential [3[ HI [5] 
is confirmed as well as the accordance with the asymptotic freedom behaviour of the Og- 
One can see that the values of model parameters rrii and are not in contradiction with 
the values of Standard Model constants such as current quark masses and the strong 
interaction coupling constant 

4. Phenomenological mass formulae for isovector gg'— mesons 

In spite of QCD difficulties which exist in for evaluations of mass spectra of hadrons 
consisting of light quarks a number of relations among hadron masses have been obtained 
using symmetry or phenomenological considerations, such as the mass relations [H] for 
Regge trajectories p[2] or the Gell-Mann-Okubo relation [131 [S]. According to the Regge 
trajectories approach a hadron with its spin J and mass M within some errors belongs 
to a straight trajectory on the (M^, J)— plane with a slope a' and a intercept cto 

J = aQ + a'M^ (9) 

Some hadrons belong to trajectories, so called, daughter trajectories, which are roughly 
parallel to the main trajectory and are distinguished with different values of a radial 
quantum number or n, n = n'^ + 1 flbi \T6\ [T7j. 

Now there is a growing interest in improving the accuracy of existing mass relations 
and obtaining these relations in the framework of the QCD or QCD inspired models. For 
instance, in Ref. [12] on the basis of the renormalization group it has been shown that 
the intercept ao cannot be calculated as a function of the coupling constant. In Ref. [IH] 
the expressions for the slope a' and the intercept has been presented in terms of the 
string tension. On the base of finite energy QCD sum rules the following mass formulae 
have been obtained for radially excited p— mesons [20] and ir— mesons |21] : 

M2(p") = M\p){2n'- + 1), (10) 
M2(7r") = M2(7r')n^ (11) 

where M'^{p) is the mass squared of the p— meson, M^(7r') is the mass squared of the first 
radial excitation of vr— meson. In the work ^22j more precise formulae for radially excited 
hadrons have been proposed. These formulae in general case can be written as 

M'^ = M',+p\n-l) (12) 

and describe trajectories on the (n,M^)— plane with different Mq values and approxi- 
mately the same p for each trajectory similarly Chew-Frautschi plots on the (M^, J)— plane 
and Nambu-Veneziano formulae for daughter Regge trajectories. 

Interesting relations for orbital and radial exited mesons' mass spectra have been 
obtained in the model of chromoelectric tubes [23] . For mesons consisting from a massive 
quark and massless antiquark the mass spectrum describes with the formula: 

= na{L + 2n + 3/2), (13) 



6 



where a is the string tension, while for mesons consisting from two massless constituents 
the mass formula is 

M"^ = 2na{L + 2n + 3/2). (14) 

Further on we restrict ourselves to those mesons which have isotopical spin values 
equal to 1, in order to bypass problems concerning unknown in some cases structure and 
mixing parameters for different mesons with the zero isotopical spin. We neglect mass 
splittings within isotopic multiplets, so systematical errors of the order of 10 MeV for 
the phenomenological scheme considered are without doubt admissible. Moreover, taking 
into account the existing experimental errors for meson masses we assume 30-j-40MeV^ as 
the typical absolute errors of the mass values evaluations, as in Section 3. 

In the framework of relativistic independent quark model a mass formula for gg'— mesons 
has a structure which differs from structures of mass formulae in other types of potential 
models. For instance, in Refs. [23] the following mass formula has been considered 

M(n2^+iLj) =i?iK,ji,c,x) + E2K,J2,c,x), (15) 

where the mass terms (mass or energy spectral functions) Ei{ni,ji,c,>i),i = 1,2, for a 
quark and an antiquark are defined as 



c + W2nr + L + j~ l/2, L + j-l/2 = 2k, k = 0,l,... 
i^[n,j,c,x) I xv/2n, + L+j-l/2, L + j - 1/2 = 2A; + 1, k = 0,l,... ^ ' 

In order to exclude superfluous meson states, the following selection rules for qq'— 
mesons with given J^^ values, quark masses mi and m2 and quantum numbers ji and j2 
are used 

ji=j2 = J+l/2, if J = L + S, 

Ji =J2 + 1 = J + 3/2, tf J^L + S, mi<m2, ^ ' 

while the radial quantum numbers n\ = = n — 1 for n'^^'^^Lj - state. The meson parity 
P = (— 1)^"*"^ and the eigenvalue C of charge conjugation for the neutral meson of qq type 
is equal to (— l)-^"*""^, where the total spin S* = or 1. 

Functions Ei{n^,ji,c, x) , i = 1,2, give the relativistic effective energies of the quark 
and antiquark moving in the mean field inside the meson, and include also an energy of 
the mean field and possible nonpotential corrections, which cannot be evaluated within 
the mean field approximation and according to the formula ([2]) should be taken into 
account with the help of the constants -Eom- Note that nonpotential corrections are the 
most important for mass spectra of light mesons [25l [26] . In general case the main part 
of mass term Ei{n1, ji,c, x), as it has been described in Section 2, are determined from 
the solution of the Dirac equation. However, for the superlight mesons, which consist of 
u— and d— quarks and antiquarks only, the phenomenological energy spectral functions 
Ei{ni,ji, c, x) in the form (fT6|) are suitable within the assumed accuracy. When evaluating 
masses of unknown excited meson states we use the formulae written above together with 
the values of two parameters c and x, which have been obtained by fitting the mass values 
of experimentally detected meson states. The c and x values obtained by this manner 
are c = 69 MeV, x = 382 ±4 MeV [27]. The results obtained with c = 69 MeV, x = 385 
MeV are shown in the Table 2, where we present the mass values of orbital excitations 
of superlight meson states up to L = 4. 
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One can obtain with the help of the formulae (ITSll . (|T6l) and the selection rules ( ITTll 
a great number of mass relations, which are fulfilled within the systematical errors of 
the phenomenological scheme considered. For instance, in the case of orbitally excitated 
vector mesons formulae f|T5|) and f|T6l) give the p— trajectory, while for radially excitated 
vector mesons they bring in the formula ffTOj) . In general case in the framework of this 
approach instead of the meson trajectories in (M^, J)— and (n, M^)— planes we have in 
(L, n, M)— space the following four series or trajectories for n^^^^Lj— meson states. 

There is a (vrp)— series for n^L^-i— meson states with P = (—1)'^+^, C = (— 1)'^+^. 
Mass values of the members of the (vrp)— series are determined by the formula: 

K'l = c + xv/2(n'^ + L) - 1 + xV2(n'^ + L) ( 18) 

There is a tt— series for n^Li— meson states with P = (— 1)^"*"^, C = (—1)'^. Mass 
values of the members of the vr— series are determined by the formula: 

= 2c + 2H^/2{n^ + L) (19) 

There is a (ptt)— series for n^L^,— meson states with P = (—1)^+-'^, C = (— l)'^^"^. 
Mass values of the members of the (pvr)— series are determined by the formula: 

M^L = c + >^^/2{n- + L) + >c^2{n^ + L) + l (20) 

There is a p— series for n^L^+i— meson states with P = (—1)-^+^, C = (—1)"^+^. Mass 
values of the members of p— series are determined by the formula: 

M^_^ = 2xV2(n'- + L) + 1 (21) 

It should be noted that only two phenomenological parameters enter in these formulae, 
which are in actual fact equal to half the pion mass and half the p— meson mass. With 
the help of these formulae for the series presented many mass relations which are not 
dependent on intrinsic quantum numbers L and n easily can be obtained. For instance, 
mass values of members of tt— series, (pn)— series and p— series with the same L obey 
the following mass relation: 

2M:iin'LL) = M:jn'LL+^) + M; Jn^Lz.) (22) 

5. Link of model of quasi-independent quarks with the model of constituent quarks 



As indicated above in order to improve the accuracy of calculations for meson mass 
spectra in the framework of quasi-independent quarks model an account of residual in- 
teractions' contributions jointly with different values for mean field energy is needed. We 
find in Sec. 3 the values of mean field energy in combination with magnitudes of spin-spin 
interaction for quarks and antiquarks with different flavours inside the pseudoscalar and 
vector mesons. Further on we propose that these values are not simply some phenomeno- 
logical constants, but represent in a rough approximation equidistant levels for mean 
field energy, which do not depend on quark flavours. We consider ordinary (non exotic) 
hadrons and relate the model of quasi-independent quarks to the constituent quark model. 
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It will be obtained that energies of fermionic constituents Ei (quark mass terms) for me- 
son ground states are equal to constituent quarks' masses within uncertainties of model 
of quasi-independent quarks. Some complications arise when one treats exited states. In 
the model of quasi-independent quarks E^ remains an exited energy of i-th constituent, 
while in the constituent quark model besides constituent quarks' masses an additional 
contribution representing energy of excitation must be taken into consideration. For the 
pseudoscalar and vector mesons the spin-spin interaction may significantly account for a 
residual interaction in the ground states of {qq') system within the accuracy of the model 
and mesons' masses can be evaluated by means of the formula ([8]) presented in Sec. 3. It 
should be noted that in the framework of the constituent quark model one can find the 
similar formula [2HI [29] : 

Mm = 1710 + 1711+1712 + (SiS2)fss, (23) 

where mi, m2 are masses of constituent quarks, mo is some additional phenomenological 
contribution. So we relate constituents' energies Ei,i = 1,2, with constituent quarks' 
masses mi,i = 1,2. In what follows this assumption will be supported by numerical fit 
for meson masses. Formulae ([8]) and fl23|) which correspond each other may be named the 
Zeldovich-Sakharov formulae. Note that constituent quarks can be applied successfully 
not only in hadron spectroscopy, but for description of effects in hadron collisions as well 
(see, e.g. [30]). 

At present due to experimental uncertainties and mixing effects there is no only one 
solution for the set of parameters entering in Eq.([8]). For instance, there is a well known 
difficulty related to the mass value of the vr— meson. Below we give an admissible fit 
with minimal number of parameters for the mass values of pseudoscalar and vector 
mesons in the ground states, which can be considered as definite [qq') systems involv- 
ing u—, d—, s—, c— , 6— quarks and antiquarks. To do this, we use in an essential manner 
the results obtained in Sec. 3. As indicated above in spite of the fact that the quark and 
antiquark energies Ei and E2 bring in a main contributions to meson masses in order to 
reduce relative uncertainties of evaluations up to 10~^ level for all mesons the additional 
contributions due to the Eq and Vgs are needed. Along similar lines Eq and Vss contri- 
butions were considered in certain of approaches [23 ESI ES] , however some ambiguities 
occur when Eq and Vgs are evaluated. When different fits for meson masses had been 
performed the case was adopted, such that the Vss value for u—, (i— quarks and antiquarks 
was equal to lOOMeV and the Eq value for p— meson was zero. In this case the Eq value 
for TT— meson is double that the £"0 value for K meson: Eq = 2Eq , where £^0 values are 
marked at the top with mesons' names. The nonzero values of Eq and Vss given in MeV 
units for different pseudoscalar and vector mesons are listed below. 

£^ = -118, E'o^/'^ = -125, EH' = -230, = -600, 
Ej = -450, El^^ = -112, = -270, * = 70, Kt = 50, (24) 
Kf = 27, = 5, ^ = 3, Vf; = 3, Kf = 2, * = 1.3, = 1 

The uncertainties of determination of the Eq and Vss values are of the order of several 
MeV. As is seen in the list of the Eq values above all nonzero Eq values for ground states 
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of pseudoscalar and vector mesons for all quark flavours are approximately multiples to 
the Eq . Thus we can accept that the Eq values are varied according to the rule: 

= -708 + 71^0^, (25) 

where n is a negative integer. To give also the values of constituents' energies Ei for 
M— , d—, s—, c—, 6— quarks and antiquarks (for two super light u— and c?— quarks and an- 
tiquarks the mean value En is presented): 

En = 337 ±3,Es = 485 ±8,Ec = 1610 ± 15, Eb = 4952 ± 20 (26) 

Notice that within the model uncertainties the constituents' energies for light quarks and 
antiquarks agree with the constituent quark masses obtained in Ref . [32] ■ Moreover if one 
take into account the difference between masses of U- and D-quarks obtained in this work, 
which is equal to 4 MeV, the masses of U- and D-quarks get the following values: 

Eu = 335 ±2, Ed = 339 ± 2 (27) 

Masses of constituent quarks fl26|) and fl27j) can be used for the derivation of such important 
parameters of Standard Model as quark mixing angles in the Fritzsch-Scadron-Delbourgo- 
Rupp approach [HU |32l [33] . 

The model of quasi-independent quarks supplemented by the proposition about the 
equidistant discrete levels of mean fleld may be named the extended model of quasi- 
independent quarks. In the framework of this model it is possible to evaluate the pseu- 
doscalar and vector mesons' masses with uncertainties of order 10~^. From the results 
presented it may be inferred that the extended model of quasi-independent quarks is a 
workable generalization of the constituent quarks model. 

6. Conclusions and discussion 

It is highly plausible that the formation of constituent quarks take place in the 
nonperturbative region with the characteristic scales, which can be evaluated with the 
universal coefficient of a slope of a linear rising potential cr = 0.20 ± O.OlGeV^^. The 
scales' values with dimensions of mass and length are equal to = 0.45 ± 0.02GeV, 
Ac = 0.44 lb 0.02Fm. As this takes place, fic defines typical magnitudes of transversal 
impulses < px > for quarks-partons inside hadrons, while a radius of a perturbative re- 
gion surrounded a current quark is equal to rc = Ac/2 = 0.22 ± O.OlFm. The region 
r > rc is most likely to be the region of a formation for a constituent quark due to non- 
perturbative interactions. If one take into account typical sizes of hadrons, then a radius 
of a constituent quark is 0.25 -v- 0.35Fm. The absolute value of Eq , which determinates 
the characteristic scale for mean field energies in the framework of the extended model of 
quasi- independent quarks can also be expressed in terms of /ic as /ic/4. 

There is a widespread opinion, that the confinement of quarks and gluons can be es- 
tablished rigorously in the QCD framework. Is this is the case, a exact dependence should 
take place, which relates Aqcd to the characteristic parameter of strong interactions in 
the nonperturbative region - the string tension a. However, in the case if the confinement 
cannot be demostrated in the QCD framework, there has to add the number of strong 



10 



interaction parameters and to introduce, at the least, one more parameter - the string 
tension. From this standpoint of some interest is investigations concerning with the gen- 
erahzation of the space-time Poincare symmetry of the QCD in the confinement region, 
for instance, to the inhomogeneous pseudounitary symmetry IU(3,1) [31], and in this case 
the universahty of the confinement potential verified in the framework of the relativistic 
model of quasi independent quarks have appresiable significance. 

In this report we considered the relativistic hadron model with quasi- independent 
quarks and presented the results of the evaluations of mass spectra of radially and or- 
bitally exited meson states. It is important to compare evaluated meson masses with 
experimentally determined ones [Tj. As it is seen from Tables 1 and 2 there is an agree- 
ment between the evaluated and experimentally determined values on the level of 1 — 2a. 
For instance, let us compare with the data p] the mass relation (l22ll . which must obey 
for the mesons with L 7^ and S = \. At present this relation can be checked only 
for the mesons with L = 1. After the substitution the experimental values of the masses 
for the ai(1260), 6i(1235) and a2(1320) mesons, we obtain, that this relation is fulfilled 
with account of experimental uncertanties. Moreover, the mass formulae presented above 
permit to explain the degeneracy with respect to mass values for different J^*" mesons. 
For example, the mass formulae ( l2Ti) give M{p") = M{p3). If one use the data, then 
1720 ± 20MeV = 1688.8 ± 2.1MeV. In the high mass value region the degeneracy of such 
kind will be increased, as it is seen in the ~ 2315MeV region, where the mass values for 
0~+, 1 , 4~+, 4 and 5 mesons are predicted. 

As it follows from the Table 2, the mass value calculated for orbital excitation of the 
vector 1 meson is equal to 1506MeV. So we expect that in the mass region between 
1.3 GeV and 1.6 GeV two standard gg'— meson resonances with J^^ = 1 and 1 = 1 
exist, namely the first radial and orbital excitations mixing each another. However, it is 
possible that the more complicated situation should be considered when in the mass region 
1.2-^ 1.8 GeV the mixing of the standard 1 mesons and the vector non gg'— mesons 
take place as well. We considered above only possible mixing between the first radial 
and orbital excitated standard 1 mesons. The most complicated situation may arise, 
which is not disscussed here, when the mixing between two radial and orbital exitations or 
with vector cryptoexotic states can occur. These cases demand further investigation. For 
instance, in the framework of the well-known potential model [10] the mass values of 2 '^5'i, 
l^Di and 3 ^5*1 states lie considerably higher (at 1.45, 1.66, 2.00 GeV, correspondingly) 
than the predictions considered. Moreover, in Ref. [35] the mass value lASGMeV has 
been evaluated the first radial excitation of p— meson. In Ref. [36] the mass values for the 
radial exitations of p— meson, which are equal to 1.4, 1.8, 2.13 Gev, have been obtained. 
Some difficulties concerning with the identification of the excited states of p— meson are 
available in Ref. [37] . 

Another importain problem is the interpretation of mesons discovered in scalar chan- 
nel. The authors of Ref.p2j found that in the region ~ IGeV there are P— wave gg'— mesons, 
The results of our evaluations (Table 2) support the existence of the P— wave gg'— mesons 
with masses ~ 980 MeV. The complementary reasoning in favour of the existence of the 
P— wave gg'— mesons with masses ~ 980 MeV is the coincidence of the evaluated mass 
value for the first radial excitation of ao(980) meson with the mass value of the ao(1450) 
meson [I]. This meson has been predicted for the first time in Ref. |3^ in the framework 
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of finite energy QCD sum rules. Although the mass values of the ao(980) and /o(980) 
mesons show the possible ideal mixing between them the problems with the intensities 
of different decay modes exist (in particular, the KK decay mode enhancement). The 
possible explanation of these facts now comprises the four-quark nature of the ao(980) 
and /o(980) mesons [391 HQ]. However, there is an accidental degeneracy of scalar mesons 
mass values with the value of i^i^— threshold, which complicates the mechanism of an 
extraction of ao(980) and /o(980) decays characteristics [H]. Moreover, it is plausible 
that the mixing of these mesons with the KK—moXecvle arises [12]. Taking into account 
the existing difficulties with identification of light mesons in the 1 -h 2 GeV mass region we 
see that allowing for the mixing between two or three mesons with different structures but 
in some cases with about the same masses is admissible for an adequate data description 
in both scalar and vector channels. 

The author is grateful to Yu.V. Gaponov, A.E. Dorokhov, A.L. Kataev, A.V. Leonidov 
V.l. Savrin, S.V. Semenov,A.M. Snigirev and V.E. Troizky for useful discussions and 
the organizers of the PNPI Winter School for the hospitality. The work was partially 
supported by the grant 26 on fundamental researches of the RRC "Kurchatov Institute" 
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Table 1. 

Evaluated masses of the 1 and mesons in comparison with the data from 
Ref. d]. 



Meson 


Mn [MeV\ 


M„ [Me 1/ J 


Meson 


Mn [Me K J 


71 /rfh r 71 /T T T"! 

M„ [Mel/J 


P 


775.5±0.4 


740 


B 


5279±0.8 


5250 


P 


1459±11 


1455 


Bs 


5367.5±1.8 


5370 


P 


1720±20 


1730 


Be 


6286±5 


6300 





1019. 460±0. 019 


1010 


Vb 


9300±40 


9330 


1 f 

<P 


1680±20 


1650 




1300±100 


1290 


1 II 






2050 


K 




1400 




3096. 916±0. Oil 


3060 


D 




2450 


1 f 


3686. 093±0. 034 


3650 


D's 




2560 


1 II 

ijj" 


4039±01 


4070 




L"\ /". c\ i~\ 1 A 

3638±4 


3600 


1 III 


^ A -1 1 ^ 

4421±4 


4390 


B' 




5650 


1 


y4oU.oU±U.ZO 


94 iu 


Bs 




/oU 


■y' 


10023.26±0.31 


9990 


B[ 




6800 


y// 


10355.2±0.5 


10325 


ib 




9960 




10579.4±1.2 


10550 


Tl" 


1812±14 


1810 




10865±8 


10830 


K" 




1850 




11019±8 


10985 


D" 




2880 


TT 


138±3.1 


120 






2950 


K 


495.65±0.02 


500 






3970 


D 


1867.7±0.4 


1850 


B" 




6070 


Ds 


1968.2±0.5 


1990 


B"s 




6130 


Vc 


2980.4±1.2 


2990 


Bl 




7150 



Table 2. 

Evaluated masses in MeV for the ground states and the orbital excitations of the qq'— mesons 
in comparison with the data from Ref.[l]. 



Meson 


jPC 




M^^(Me1/) 


Meson 


JPC 


M^^P(MeV")^i^l 


M^^(MeV") 


TT 


0-+ 


138.039±0.004 


138±30 


P3 


3— 


1688.8±2.1 


1722±30 


P 


1"- 


775.5±0.4 


770±30 


al 


2++ 




1873±30 


ao 


0++ 


984.7±1.2 


998±30 


bl 


3+- 




2024±30 


hi 


1+- 


1229.5±3.2 


1227±30 


al 


3++ 




2031±30 


ai 


1++ 


1230±40 


1280±30 


a^ 


4++ 


2001±10 


2037±30 




2++ 


1318.3±0.6 


1334±30 


al 


3— 




2177±30 


a\ 


1 — 




1506±30 




4-+ 




2316±30 




2-+ 


1672.4±3.2 


1678±30 


aj 


4— 




2313±30 


al 


2— 




1700±30 


4 


5— 




2310±30 
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